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CENTRAL LIMIT THEOREM FOR FOURIER TRANSFORMS OF 
STATIONARY PROCESSES 

By Magda Peligrad 1 and Wei Biao Wu 2 

University of Cincinnati and University of Chicago 

We consider asymptotic behavior of Fourier transforms of sta- 
tionary ergodic sequences with finite second moments. We establish 
a central limit theorem (CLT) for almost all frequencies and also an 
annealed CLT. The theorems hold for all regular sequences. Our re- 
sults shed new light on the foundation of spectral analysis and on the 
asymptotic distribution of periodogram, and it provides a nice blend 
of harmonic analysis, theory of stationary processes and theory of 
martingales. 

1. Introduction. In frequency or spectral domain analysis of time series, 
periodograms play a fundamental role. Since its introduction by Schuster 
(1898), periodograms have been used in almost all scientific fields. Given 
a realization (X,)™ =1 of a stochastic process (Xj)j^z, the periodogram is 
defined as 

1 n 

^Xjexpijit 

3=1 



G 



where i = y/—l is the imaginary unit. Periodogram is the building block in 
spectral domain analysis and a distributional theory is clearly needed in the 
related statistical inference. If (Xj) is a Gaussian process, then the Fourier 
transform 

n 

S n (9) = ^2x j exp(ji9) 

3=1 
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is complex Gaussian. Fisher (1929) proposed a test for hidden periodicities 
and obtained a distributional theory based on i.i.d. Gaussian random vari- 
ables. If (Xj) is not Gaussian, the distribution of S n (9) typically does not 
have a close form and one needs to resort to asymptotics. It is well known 
since Wiener and Wintner (1941) [see also Lacey and Terwilleger (2008)] 
that for any stationary sequence (Xj)j G z hi C 1 (namely i^Xo] < oo) there 
is a set f2' of probability 1 such that for all 6 and uj S 0', S n (9)/n converges. 
Our problem is to investigate the speed of this convergence by providing a 
central limit theorem for the real and imaginary parts of S n {6) / sfn. 

The above central limit problem was considered by many authors under 
various dependence conditions. We mention Rosenblatt [(1985), Theorem 
5.3, page 131] who considered mixing processes; Brockwell and Davis [(1991), 
Theorem 10.3.2, page 347], Walker (1965) and Terrin and Hurvich (1994) 
discussed linear processes, and Wu (2005) treated mixingales. Other con- 
tributions can be found in Olshen (1967), Rootzen (1976), Yajima (1989), 
Woodroofe (1992), Walker (2000), Lahiri (2003) and Lin and Liu (2009) 
among others. 

To establish an asymptotic theory for S n (9)/y/n, we shall provide the 
framework of stationary processes that can be introduced in several equiv- 
alent ways. We assume that (£ n )nez is a stationary ergodic Markov chain 
defined on a probability space (£l,J-,P) with values in a measurable space. 
The marginal distribution is denoted by ir(A) = -P(£o £ A). Next let Cq(tt) 
be the set of functions such that J h? dir < oo and J hdir = 0. Denote by 
•7-fc the ex-field generated by £j with j < k, Xj = h{^j). For any integrable 
random variable X we denote E^{X) = E^XlJ 7 ^). We assume h € Cq(it); in 
other words we assume ||Xo|| := (E\Xq j) 1 / 2 < oo and E(Xq) = 0. Notice that 
any stationary sequence (Yk)k^i can be viewed as a function of a Markov 
process = (Yj-,j < k) with the function g(£,k) = ^k- 

The stationary stochastic processes may be also introduced in the fol- 
lowing alternative way. Let T:f2 1— > f2 be a bijective bi-measurable trans- 
formation preserving the probability. Let J-q be a cr-algebra of J- satisfy- 
ing Fx C T~ 1 (J r o). We then define the nondecreasing filtration {F^j^z by 
JFj = T~ J (J : o) (referred to as the stationary filtration). Let Xq be a random 
variable which is J^-measurable. We also define the stationary sequence 
(Xj)j^z by Xj = Xq o T 3 . In this paper we shall use both frameworks. 

The rest of the paper is structured as follows. The main results are pre- 
sented in Section 2 and proved in Section 4. Our proofs in Section 4 provide 
an interesting blend of harmonic analysis, martingale approximation and 
theory of stationary processes. Examples of regular processes and further 
extensions are given in Section 3. 
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2. Main results. We shall assume that the following regularity condition 
holds: 

(2.1) ^(Xol-F-oo) = 0, P-almost surely, 

and also that the sequence is stationary and ergodic. The regularity con- 
dition is quite mild and it is satisfied for many popular processes used in 
practice. Section 3 provides examples of stationary ergodic processes for 
which (2.1) holds. 

We shall present first a central limit theorem for almost all frequencies. 
In Theorem 2.1, we let the parameter 6 be in the space [0, 2tt], endowed 
with Borelian sigma algebra and Lebesgue measure A. We denote by "=>•" 
the weak convergence, or convergence in distribution. 

Theorem 2.1. Let (Xk)k<=z be a stationary ergodic process such that 

(2.1) is satisfied. Then for almost all 9 £ (0, 2ir), the following convergence 
holds: 

(2.2) hm E \ SnW = g {9) (say), 

n— >oo n 

where g is integrable over £ [0,2-7r], and 

(2.3) ^=\Re(S n (0)),1m(S n (0))] => [N 1 (e),N 2 (9)} under P, 
\/n 



where N\{0) and Ni{6) are independent identically distributed normal ran- 
dom variables mean and variance g(9)/2. 

As implied by Lemma 4.2 in Section 4, g(9)/(2ir) is actually the spectral 
density associated with the spectral distribution function induced by the 
covariances 

(2.4) c,=cov(X ,X i ), jeZ. 

More specifically, by Herglotz's theorem [Brockwell and Davis (1991)], there 
exists a nondecreasing function G (the spectral distribution function) on 
[0, 2tt] such that, for all j £ Z, 

/•27T 

(2.5) Cj = / exp(ij9)dG(9). 

Jo 

Hence, by Lemma 4.2, G is absolutely continuous and the spectral density 
G'{9) equals to g{6)/{2ir) almost surely. By (2.5) or (4.3), f Q 2n g(9) d9 = 2irc Q . 
So a nice implication of our results is that, under the regularity condition, 
we obtain an interesting representation of the spectral densities (see Lemma 
4.2 for details). 
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Following the proof of Theorem 2.1, by the Cramer- Wold device, for 

\ n 



we have that, for almost all pairs (9', 6") (Lebesgue), V n (uj,6') and V n (uj,9") 
are asymptotically independent. In this sense Theorem 2.1 justifies the folk- 
lore in the spectral domain analysis of time series: the Fourier transforms 
of stationary processes are asymptotically independent Gaussian. Namely, 
in the spectral or Fourier domain, the Fourier-transformed processes are 
asymptotically independent, while the original process can be very strongly 
dependent (see Example 3.3). 

Theorem 2.1 substantially improves the result in Wu (2005) that proves 
(2.3) under the following stronger condition: 



(2-6) E 



\\E{X n \F»)f 



< oo. 



n 

71=1 



We shall also establish the following characterization of the annealed CLT. 
Let Id2 denote the identity 2x2 matrix. 

Theorem 2.2. Under the same conditions as in Theorem 2.1 on the 
product space ([0, 2tt] x x J, A x P) we have 

(2.7) -^[Re(S n (9)),lm(S n (9))]^[g(U)/2] 1 / 2 N(0,ld2) under X x P. 



n 

Here U is a random variable independent of iV(0, Id2) and uniformly dis- 
tributed on [0,27r] and g{-) is defined by (2.2). 

Two types of stochastic processes can be considered concerning the partial 
sum S n (9). The process V n (u,9) indexed by 9 is asymptotically Gaussian 
white noise. For another version, we consider 

S 

W n (t,Gj,9) = -±p-(u),9), 0<t<l, 



n 



where [^J = max{fc 6Z:fc<i} is the integer part of x. We shall prove the 
following invariance principle: 

Proposition 2.1. Assume that (X^) is stationary ergodic and satisfies 
(2.1). Then W n (t,u,9) is tight in 13(0,1), and 

[Re(W n (t,cj,9)),lm(W n (t,cj,9))} 

=► [g(U)/2}^ 2 [W'(t),W"(t)} under A x P, 

where (W'(t),W"(t)) are two independent standard Brownian motions inde- 
pendent of U , and U is a random variable uniformly distributed on [0, 2n] . 
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Remark 2.1 (Nonadapted case). Our CLT also holds if Xq is not JFq- 
measurable. Then clearly the regularity condition we shall impose is 

Xq is J^-measurable and E{Xq\J : -. O0 ) = almost surely. 

Remark 2 . 2 (Adapted nonregular case) . For general adapted sequences 
our CLT result still holds under centering. Let S n (9) = S n {6) - E(S n (6)\F- 00 ) 
and Xk = Xk — E{X k \J : - 00 ), where E(Xi : \J r - 00 ) denotes the following limit 
that holds almost surely and in C 2 

lim E(X k \T. n ) = E(X k \T^ 00 ). 

n— voo 

Then, for almost all £ (0,27r), we have 

-^=[Re(S n (0)),Im(S n (0))} ^]W0,^Id 2 J under P 



by applying Theorem 2.1 to the stationary sequence X k . Therefore the con- 
clusion of Theorem 2.1 holds if we replace the assumption of regularity (2.1) 
by the following: for A-almost all 8 

—=E{S n {9)\F^ oa ) — > in probability. 



Now, since \\E(S n {e)\F-oo)h < \\E(S n {e)\F^ n )\\ 2 < \\E(S n {6)\? Q )\\ it Theo- 
rem 2.1 still holds if 

" \\E(S n {e)\T- n )\\ 2 ^Q asn^oo 
n 



or under the condition 

(2.8) -±=\\E{S n {e)\Fo)h^Q asn^oo. 

Jn 



Remark 2.3 (Conditional CLT). Since we use in the proof martingale 
approximation actually our CLT is a conditional CLT, that allows for a 
random change of measure. See Hall and Heyde (1980) and Dedecker and 
Merlevede (2002). 

Remark 2.4 (Periodogram). We notice that, as a consequence of Theo- 
rem 2.1, for sequences satisfying (2.1) the periodogram n~ l \S n {6)\ 2 is asymp- 
totically distributed as ^4^x 2 (2) for almost all frequencies. 
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Remark 2.5 (Resulting identities). By using the invariance principle in 
Proposition 2.1 we can get the convergence of many interesting functionals 
of S n (6) and periodograms. As a consequence of Proposition 2.1 we can get, 
for instance, 



1 



|X || 2 n7r J 
-^E 



E 



de 



max > Xucos(k9) 
- - fc=l 

sup W{t)\ = / 2[l-$(Vy)]dy = l 

<t<l 1 JO 



by noting that P(sup < t < 1 W(t) >u) = 2P(W(1) > u) for u > 0. Here $( 
is the standard Gaussian distribution function. 



3. Examples. Here we present several examples of processes for which 
the conclusions of Theorems 2.1 and 2.2 hold. 

Clearly condition (2.1) is satisfied if the left tail sigma field P-oo is triv- 
ial. These processes are called regular [see Chapter 2, Volume 1 in Bradley 
(2007)]. Notice, however, that our condition (2.1) refers rather to the func- 
tion Xq = /(£o) m relation to the tail field J 7 ^. 

Example 3.1 (Mixing sequences). We shall introduce the following mix- 
ing coefficients: for any two a-algebras A and B define the strong mixing 
coefficient 

a{A, B) =sup{|P(^ n B) - P{A)P(B)\ :AeA,BeB} 

and the p-mixing coefficient, also known as maximal coefficient of correlation 

p(A,B) = sup{Cov(X,Y)/\\X\\ 2 \\Y\\ 2 : X € C 2 (A),Y G C 2 (B)}. 

For the stationary sequence of random variables (X^kez, F n denotes the 
cr-field generated by Xi with indices i>n, and T m denotes the cr-field gen- 
erated by Xi with indices i <m. The sequences of coefficients a(n) and p(n) 
are then defined by 

a(n) = a(Po,P n ) and p{n) = p(P ,P' n '), 

respectively. For strongly mixing sequences, namely the strong mixing coef- 
ficients a(n) — > 0, the tail sigma field is trivial [see Claim 2.17a in Bradley 
(2007)]. Examples of this type include Harris recurrent Markov chains. If 
lim n ^ 00 p(n) < 1, then the tail sigma field is also trivial [see Proposition 5.6 
in Bradley (2007)]. 

Example 3.2 (Functions of Gaussian processes). Assume (Y^) is a sta- 
tionary Gaussian sequence and define X n = /(Y&, k < n). Let / be such that 
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E(Xq) = and E(Xq) < oo. Since any Gaussian sequence can be repre- 
sented as a function of i.i.d. random variables, the process is then regular. 
Rosenblatt (1981) considered Fourier transforms of functionals of Gaussian 
sequences. 



Example 3.3 (Functions of i.i.d. random variables). Let be i.i.d. and 
consider X n = f(sk, k <n). These are regular processes and therefore The- 
orems 2.1 and 2.2 are applicable. Examples include linear processes, func- 
tions of linear processes and iterated random functions [Wu and Woodroofe 
(2000), among others]. For example, let X n = Yl'jLo a j 
with mean and variance 1, and aj are real coefficients with 
In this case X n is well defined, and, by Lemmas 4.1 and 4.2, the spectral 
density is g(6)/(2ir), where 



j, where Ej are i.i.d. 
'Vf <oc. 



9(0) 



As a specific example, let aj = j ' /logj, j > 2, and ao = ai = 1. By el- 
ementary manipulations, the covariances Cj ~ (logj)^ 1 , which decays very 
slowly as j — > oo, hence suggesting strong dependence. For this example, 
condition (2.6) is violated. By the Tauberian theorem, as 9 — > 0, g(9) ~ 
7r/(|0| log 2 |#|), which has a pole at = 0. 

Example 3.4 (Reversible Markov chains). As before, let £j be a sta- 
tionary ergodic Markov chain with values in a measurable space. We use the 
notation and constructions from the Introduction. The marginal distribution 
and the transition metric are denoted by tt(A) = -P(£o £ A) and Q(£,o,A) = 
P(£i 6 ^4|£o)- In addition Q denotes the operator Qf(£,) = f f(z)Q(^,dz). 
Let Q* be the adjoint operator of the restriction of Q to C 2 (ft) and assume 
Q = Q* . Then, for any / £ £q(-7t) the central limit theorem of Theorem 2.1 
holds. To see this we shall verify condition (2.8). By spectral calculus 



\\E(s n (e)\?o)\? 



£(texp(tf))' 



fe=i 



Pfidt), 



where pf denotes the spectral measure of / with respect to Q [see, e.g., 
Borodin and Ibragimov (1994) for this identity]. For 9 ^ 0, 7r and — 1 < t < 1, 
we have 

2 



k=l 



<4|l-exp(i6»)t|- 2 

= 4(1 + t 2 - 2tcos9)~ 1 < 4(1 - (costf) 2 )^ 1 . 
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Therefore, for A-almost all 9 



1 



-\\E(S n (9)\T )\\^0. 
n 

4. Proofs. We shall establish first some preparatory lemmas. The almost 
sure convergence in Lemma 4.1 was shown by Wu (2005). The convergence 
in C 2 is new here. For k G Z we define the projection operator by 



(4.1) 



P k -=E(-\T k )-E(-\T k -i). 



Lemma 4.1. Let 

n 

T n (9) = Xj exp(j^) = X + S n (9). 

3=0 

Under (2.1), for X-almost all 9 (Lebesgue), we have 

oo 

(4.2) V T n {9)^^2v X l ex.p{U9)=:D {9), P-almost surely and in C 2 . 

1=0 

PROOF. By (2.1), J2 k& \\V k X \\ 2 = \\X \\ 2 2 < oo, we have £ feeZ \V X k \ 2 < 
oo, P-almost surely. Therefore by Carleson's (1966) theorem, for almost 
all u, ^2i <k<n ('PoX k ) ex.p(ik9) converges A-almost surely, where A is the 
Lebesgue measure on [0,2-/r]. Denote the limit by Dq = Dq(9). We now con- 
sider the set 

A = {(9,u) C [0, 2tt] x J], where {T'oS n (9)} n does not converge} 

and notice that almost all sections for u fixed have Lebesgue measure 0. 
So by Fubini's theorem the set A has measure in the product space and 
therefore, again by Fubini's theorem, almost all sections for 9 fixed have 
probability 0. It follows that for almost all 9, Vo(S n (9)) — > Dq almost surely 
under P. Next, by the maximal inequality in Hunt and Young (1974), there 
is a constant C such that 

\ S ^\V {S n {9))\ 2 ]\{d9)<Cy^\V Q X k \l 
1 n J k 



o 

and then we integrate 

-2-n 

E 

'o 

Therefore, 



E 



sup|P o (Sn(0))| 2 KdO) < C\\X \\l < oo. 

n 

sup|Po(<Sn(#))| 2 < oo for almost all 9. 
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Since \VoS n (0)\ < sup n \VoS n (0)\, and the last one is integrable for almost all 
9, by the Lebesgue dominated convergence we have that Vo{S n {9)) converges 
in C 2 . □ 

Lemma 4.2. Let g(9) = E\D Q (9)\ 2 . For all j G Z, we have 

<-2- 



(4.3) 



g{9) exp(jiO) dO = 2irc 



o 



3i 



where Cj are defined by (2.4)- So (cj) are the Fourier coefficients of g. Ad- 
ditionally, for almost all 9, 



(4.4) 



lin, *i!! = 9 ( 9 ). 



n 



Proof. Without loss of generality we let j > and j < n. As before, let 

n 

V Q T n {p) = Y^V Q X l exp{UO). 



1=0 

By elementary trigonometric identities, we have 

r 2n n ~3 



(4.5) 



1 

2^ 



o 



\V T n (9)\ 2 ex V (ji9)d9 = Y^PqXi){VqX u 



1=0 



Since Xj = ^Zi^ViXj, by orthogonality of martingale differences and sta- 
tionarity we have that 



Cj = lim E 

J TV-s-oo 



E w> E 

=-N / \l=-N 

lim ]T 







N-too 



l=-N 
n-j 



lim 



Z=0 



By (4.5) and the Lebesgue dominated convergence theorem, as in the proof 
of Lemma 4.1, (4.3) follows in view of Hunt's maximal inequality since 
supJ'Po x T n {9)\ is integrable. 

Now we prove (4.4). By stationarity, we have 

^ 1 n n 

-E\S n (9)\ 2 = - E E E ( X i X i) expM) exp(-iZ0) 
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1 n— 1 n 

(4.6) =-^^ Cj _ z exp((j-/)^) 



n 

3=1 1=1 
n-l 

j=-(n-l) 



— )cjexp(ij6). 



Namely E\S n {6)\ 2 /n is the Cesaro average of the sum Y^j=-i c j exp(ijf). 
Note that 5(6*) = ||L» (6 I )|| 2 is integrable over [0, 27r]. Therefore by the Fejer- 
Lebesgue theorem [cf. Bary (1964), page 139 or Theorem 15.7 in Champeney 
(1989)], (4.4) holds for A-almost all Q G [0,2tt] (Lebesgue). □ 

Remark 4.1. In the proof of Lemma 4.2, (4.6) implies that the se- 
quence (cj exp(ij#)) is Cesaro summable. It turns out that, generally speak- 
ing, Xljt=o c i x ex P(*i^) m ay not exist for almost all 9. Consider the example 
in Kolmogorov (1923) [see also Theorem 3.1, page 305, in Zygmund (2002)]: 
there exists a sequence of nonnegative trigonometric polynomials f n with 
constant term 1/2, a sequence of positive integers qu — > 00 and a positive 
sequence A n — > 00 , such that the function 



k=l ^rife 

is integrable. However, for almost all 9, the Fourier sum c i exp(/i#) di- 
verges, where q = J Q 27r g{9) exp(U9) d6 is the Fourier coefficient of g. Let 
G{x) = j Q g{u) du. By Herglotz's theorem [Brockwell and Davis (1991)], 
there exists a stationary process (X,) such that its spectral distribution 
function is G and its covariance function is q. 

Lemma 4.3. Assume (2.1). On the product space ([0, 27r] x S~l,B x J 7 , A x 
P) we have that 

'max 1<k<n \S k (9)\ 2 \ . 

is uniformly integrable. 



n >n>l 



Proof. Let m be a positive integer. We shall decompose the partial 
sums in a sum of m martingales and a remainder in the following way: 

Notice that for any < j < m — 1, 
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is a martingale adapted to the filtration (B x Fk). Moreover, since (Xk)kez 
is a stationary sequence with variables square integrable, it follows that 
(X|)fc £ z is a uniformly integrable sequence. This fact implies that for j 
fixed the sequence (Vk-j(Xk))kez is also uniformly integrable. It follows 
that Yli=i exp(i/0) x Vi-j(Xi) is a martingale with uniformly integrable dif- 
ferences under the measure A x P. It is known that for a martingale with 
uniformly integrable differences we have 

k 2 

£exp(»Z0)7Vi(*i) 



l 

— max 

n Kk<n 



is uniformly integrable [see, e.g., Dedecker and Rio (2000), Proposition 1]. 
The result follows since by Hunt and Young (1974) maximal inequality 



2?r 



max 

Kk<n 



i=i 



d6<cY)Ei-m{Xi 



and therefore, denoting by E the expected value with respect to A x P, we 
have by regularity condition (2.1) that 



1 



E 



n 



max 

Kk<n 



J2^M^)El-m(Xl) 



1=1 



< — E\Ei_ m (Xi 



i=i 



CE\E_ m (X ) 







as m — > oo uniformly in n. So, the uniform integrability follows. □ 

4.1. Proof of Theorem 2.1. The first assertion of Theorem 2.1 is just 
Lemma 4.1. We now prove (2.3). 

Step 1. The construction of martingale. 

Define the projector operator by (4.1). Then we construct as in Lemma 
4.1 

n 

V l (S n {e)) = E{S n (e)\F l )-E{S n {e)\T Q )=Y j eMike)'Pi(Xk) 

k=l 

and then by Lemma 4.1 for almost all 6 

Pi(S n (e))-^D 1 (e) in£ 2 . 
To verify it is a martingale we start from 

E(Vi(S n (6))\J r o) = almost surely under P 
and by the contractive property of the conditional expectation 

= E(Vi(S n (e))\F ) -> E(D 1 (9)\T ) in c 2 . 
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We then construct the sequence of stationary martingale differences {Dk(0))k>i, 
given by 



V k {S n+k (6) - S k {9)) -> exp(ik9)D k (9) in C 2 . 

Step 2. Martingale approximation. 
Denote by 

M n {9) = ]T exp(ik9)D k (9). 

\<k<n 

We show that, for almost all 9, 

(47) E\S n (9)-M n (9)\ 2 ;a 

n 

To this end, note that S n (9) — E(S n (9)\ F$) and E{S n {9)\Fo) are orthogonal, 
we have 

(4.8) \\S n {9)\\ 2 = \\Sn{9) - E(S n (9)\F )f + \\E(S n (9)\F )\\ 2 . 

For those 9 such that (4.2) holds, we have, by the orthogonality of martingale 
differences and the stationarity, that 

\\S n (9) - E(S n (9)\F ) - M n (9)\\ 2 

n n 

(4.9) = ^2\\V k (S n (9) - M n {9))\\ 2 = \\V k S n {9) - e* ke D k {9)\\ 2 

k=l k=l 

n 

= J2W V oT n -k(0) ~ D (9)\\ 2 =o(n). 

k=l 

Hence, by (4.8) and (4.4), we have 

\\E{S n {9)\m\ 2 y ||5 n (g)|| 2 -||M n .(g)|| 2 
hm sup = hm sup = 

n— »oo Tl n—too Tl 

by noting that \\M n {9)\\ 2 = n\\D (9)\\ 2 . Hence we have (4.7) in view of (4.9). 
Step 3. The CLT for the approximating martingale. 

It remains just to prove central limit theorem for complex valued martin- 
gale 



)= exp(ik9)D k (9). 



n 

Kk<n 



As a matter of fact we shall provide a central limit theorem for the real part 
and imaginary part and show that in the limit they are independent. The 
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proof was carefully written down in Wu (2005). By the Cramer- Wold device 
we have to study the limiting distribution of the martingale 

-= V [sRe(eMik9)D k (9)) + tlm(exp(ik9)D k (9))]. 



l<fe<n 

By the Raikov-type of argument, in order to prove the CLT we have only to 
show 

- ^2 [sRe(exp{ik6)D k (e)) + t Im{exp(ik9)D k (9))} 2 -> v ^ + , 

l<k<n 

This follows from combining the following two facts. First by stationarity 
- V \exp(ik9)D k (9)\ 2 = - V \D k (9)\ 2 ^ E\D (9)\ 2 

fl — * n — ^ 



n 



l<k<n l<k<n 

and then by Lemma 5 in Wu (2005) for almost all 9 

i ]T [exp(ik9)D k (9)} 2 ^0. 



l<k<n 

The rest is simple algebra. 

4.2. Proof of Theorem 2.2. Consider the product space ([0, 2ir] xfl,B x 
J-,\ x P). Let P = A x P and E the corresponding expected value. We 
already have shown in Theorem 2.1, that for A-almost all 6, 



Eexp 



n 



(sRe(S n (9))+tlm(S n (9))) 



exp 



(s 2 + t 2 )g(6) 



for s,t € K. Then we integrate with 9 and by the dominated convergence 
theorem we obtain 



Eexp 



- 7 =(sRe(S n (9)) +tlm(S n (9))) 



1 f 2n 

Tvl exp 



(s 2 +t 2 )g(9) 



d9. 



We then identify the limiting distribution as being a mixture of two indepen- 
dent random variables: a standard normal variable with a variable uniformly 
distributed on [0, 27r]. 

4.3. Proof of Proposition 2.1. It is easy to see that the finite-dimensional 
distributions are convergent. So we just have to prove tightness. From Billings- 
ley (1999), stationarity and standard considerations this follows by Lemma 
4.3. 
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